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ABSTRACT: (Type single-spaced below line. Lettered paragraphs (a) through (d) 
should include: a. brief statement of the overall objectives and justification of the work; b. 
brief statement of the accomplishments of the prior year, or "new proposal;" c. brief listing 
of what will be done this year, as well as how and why; and d. one or two of your recent 
publications relevant to the proposed work.) 


(a) The broad objective of this work is to improve understanding of the internal 
structures and thermal and stress histories of the solid planets by combining results 
from analytical and computational modeling, and geophysical data analysis of 
gravity, topography and tectonic surface structures. 


(b) During the past year we performed two quite independent studies in the attempt 
to explain the Mariner 10 magnetic observations of Mercury. In the first we revisited 
the possibility of crustal remanence by studying the conditions under which one 
could “break symmetry” inherent in Runcorn’s model of a uniformly magnetized shell 
to produce a remanent signal with a dipolar form. In the second we applied a thin 
shell dynamo model to evaluate the range of intensity/structure for which such a 
planetary configuration can produce a dipole field consistent with Mariner 1 0 results. 

(c) In the next full proposal cycle we will: (1) develop numerical and analytical and 
models of thin shell dynamos to address the possible nature of Mercury’s present- 
day magnetic field and the demise of Mars’ magnetic field; (2) study the effect of 
degree-1 mantle convection on a core dynamo as relevant to the early magnetic field 
of Mars; (3) develop models of how the deep mantles of terrestrial planets are 
perturbed by large impacts and address the consequences for mantle evolution; (4) 
study the structure, compensation, state of stress, and viscous relaxation of lunar 
basins, and address implications for the Moon’s state of stress and thermal history 
by modeling and gravity/topography analysis; and (4) Use a three-dimensional 
viscous relaxation model for a planet with generalized vertical viscosity distribution 
to study the degree-two components of the Moon’s topography and gravity fields to 
constrain the primordial stress state and spatial heterogeneity of the crust and 
mantle. 

(d) Papers of particular relevance to the proposed investigation: 

Stanley, S., J. Bloxham, W.E. Hutchison, and M.T. Zuber, Thin shell dynamo models 
consistent with Mercury’s weak magnetic field, Earth Planet. Sci. Lett., in press, 2005. 
Aharonson, O., M.T. Zuber and S.C. Solomon, Crustal remanence as a source for 
Mercury’s magnetic field, Earth Planet. Sci. Lett., 6945, doi: 10.1016/ 
j.epsl. 2003. 11.020, 2004. 



2 


I. INTRODUCTION 

This report summarizes progress in our 
research effort in planetary geophysics 
and tectonics. During the past year our 
research group has addressed a range of 
questions that involve aspects of these 
processes as applied to the terrestrial 
planets, with an emphasis on the Moon 
and Mercury. 

During the past year we published or 
submitted or contributed to six manuscripts 
and one book chapter under the auspices 
this grant. A list of these publications 
follows the technical report. A discussion of 
future work was presented in the full 
proposal to the PGG program. 

II. PROGRESS TOWARDS 
UNDERSTANDING THE MAGNETIC 
SIGNATURE OF MERCURY 

Some Background 

Mercury represents a key towards 
understanding the evolution of the terrestrial 
planets. In striking contrast to the Moon's 
depletion in iron and small (if any) core, 
Mercury's size and mass [Anderson et al., 
1987; Anderson et al., 1996] indicates a 
high metal/silica ratio and a metallic mass 
fraction of about twice that of the Earth, 
Venus and Mars. The uncompressed 
density (5500 kg m^) suggests that if the 
planet differentiated into a silicate mantle 
and iron core then ~0.75 R^, and the 
fractional core mass is about 0.65 [Siegfried 
and Solomon, 1974]. This unusual internal 
structure combined with the puzzling 
detection of a dipole magnetic signature 
[Connerney and Ness, 1 988] during two 
Mariner 1 0 flybys has led to considerable 
debate about internal structure, particularly 
core state. The possibility that Mercury's 
magnetic field is a consequence of a core 
dynamo would require that the planet's 
metallic core is at present at least partially 
molten. For thermal evolution models in 
which core-mantle differentiation occurred 
early and the core is either pure iron or an 
iron-nickel solid solution, an initially molten 
core should have frozen out by now 
[C as sen et al., 1976; Fricker et al., 1976; 
Siegfried and Solomon, 1974], However, a 
present-day fluid core is possible if a lighter 
element such as sulfur is present [Schubert 
etal., 1988], Alternatively, perhaps the 
magnetic signature detected by Mariner 10 


indicates a frozen-in crustal field associated 
with thermal remanence, which could be 
consistent with a solid or non-convecting 
core. The evolution of core state with time 
has implications for orbital evolution: the 
presence of a fluid core at the time that 
Mercury entered its 3:2 spin orbit resonance 
would have enhanced its capture 
probability [Peale, 1988], 

Mariner 10 Observations 

Mariner 1 0 observed Mercury’s 
magnetic field during 2 flybys of the planet 
in 1974-1975 (Figs. 1 and 2), revealing the 
presence of a magnetic field with dipole 
moment of about 300 nT-R M 3 (1 R M = 2440 
km) [Ness etal., 1975; Ness etal., 1976]. 
The magnetic observations can be 
consistent with a present-day dynamo 
[Connerney and Ness, 1988], but other 
interpretations of the data are possible. 
Crustal remanence is possible but was 
discounted for some time due to a well 
known theorm indicating that a uniformly 
magnetized shell in the presence of an 
internal source will have no external field 
subsequent to the removal of the source 
[Runcorn, 1975a; Runcorn, 1975b]. 



Determining unambiguously whether 
the observed field is due to crustal 
remanence, an active dynamo, or 
thermoelectric currents is difficult 
[Aharonson etal., 2004; Giampieri and 
Balogh, 2002; Schubert et al., 1988; 
Stanley et al., 2005; Stevenson, 1987] 
because of the field’s magnitude and the 
limited spatial and temporal resolution of 
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the current data [Connerney and Ness, 
1988; Ness, 1979], For example in the last 
funding cycle we developed an analytical 
theory [Aharonson etal., 2004] (Appendix 
2} that showed that a remanent signature 
on Mercury can have a significant dipole 
component due to the latitudinal influence 
of surficial heating on the depth to the 
Curie isotherm. 

Recent ground-based observations of 
Mercury's forced iibrations in longitude 
provide compelling indirect evidence that 
Mercury’s core is at least partially fluid 
[Margot et at., 2004] ; hence a basic 
necessary condition for dynamo action 
appears to be fulfilled. However energetic 
and magnetostrophic balance arguments 
[Schubert et al., 1988; Stevenson, 1987] 
show that a dynamo source for Mercury’s 
observed magnetic field is problematic if 
one expects an Earth-like partitioning of 
toroidal and poloidal components of the 
field. 

The Future 

Future observations from the NASA 
MESSENGER mission [Solomon etal., 
2001] will provide a range of geophysical, 
geochemical and geological observations 
relevant to addressing the nature of 
Mercury's thermal evolution, with detection 
of the core state and the mechanism of 
magnetic field generation being high priority 
science objectives. In the mean time, we 
develop a suite of models to potentially 
explain Mercury’s magnetic signature. 

III. MODELS TO EXPLAIN THE 
MARINER 10 MAGNETIC SIGNATURE 
OF MERCURY 

A Remanent Magnetization Model 

Previous attempts to explain Mercury’s 
magnetic field as a consequence of 
remanent magnetization were dismissed 
[Stephenson, 1976] because of an 
assertion of Runcorn [Runcorn, 1975a; 
Runcorn, 1975b], that lacking any lateral 
variations in shell thickness an external 
magnetic field vanishes. But if the 
symmetry of the shell can be broken, then 
remanent magnetization should be 
possible. 

In the past year we investigated how 
variations in the thickness of a surficial layer 
that is available to be magnetized might be 
responsible for external magnetic fields. 

Our work [Aharonson et al., 2004] provides 


a general solution to the variable layer- 
thickness problem, demonstrates some 
special cases that are easily obtained from 
it, and applies the formulation to Mercury. 
Our aim was not to dispute that Mercury's 
magnetic field may indeed originate in the 
core, but rather to reexamine the often 
dismissed [Stephenson, 1 976] that it 
originates in the crust. 

We considered the magnetic field of a 
shell uniformly magnetized by an internal 
dipole that is subsequently removed. The 
Gauss coefficients of the resulting field were 
given in terms of the spherical harmonic 
coefficients of the shell thickness. This 
general solution can easily be reduced to 
common special cases by superposition. 
For a shell of constant thickness the 
external field vanishes (by Runcorn's 
theorem). But for a laterally varying 
temperature field, such as would be 
expected for Mercury due to latitudinal 
differences in illumination and longitudinal 
differences associated with Mercury’s orbit, 
the resulting magnetic moments are 
appreciably greater than the previously 
published correction due to rotational 
flattening. We showed that if the crust of 
Mercury contains rocks capable of 
sustaining high specific magnetizations, 
then the Mariner 10 observations of 
Mercury's magnetic field are consistent in 
magnitude and geometry with the 
predictions of this model [Aharonson etal., 
2004], For such a scenario, the requirement 
of a fractionally large molten outer core 
would be relaxed. . 


A Thin Shell Dynamo Model 

In a preliminary study, we [Stanley et 
al., 2005] used a formulation for a 3-D 
numerical dynamo model [Kuang and 
Bloxham, 1997; Kuang and Bloxham, 

1999] to demonstrate that if Mercury’s core 
consists of a thin fluid shell surrounding a 
solid core (the geometry suggested by 
some thermal evolution models for Mercury 
[Schubert et al., 1988; Stevenson, 1987; 
Stevenson etal., 1983], then a thin shell 
dynamo is capable of producing fields with 
toroidal-poloidal field partitioning similar to 
Mercury (and different from Earth). The 
purpose of the study was to determine 
whether dynamo models capable of 
explaining Mercury’s observed magnetic 
field plausibly could have existed. 
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As shown in Fig. 3, we examined the 
ratio of the dipole field at the core-mantle 
boundary to the toroidal field in the core for 
various shell thicknesses and Rayleigh 
numbers. We found that some thin shell 
dynamos can produce magnetic fields with 
Mercury-like dipolar field intensities. In 
such dynamos, the toroidal field is 
produced more efficiently through 
differential rotation than the poloidal fhe 
poloidal field is produced through 
upwellings interacting with the toroisal field. 
The poloidal field is also dominarted by 
smaller-scale structure that was not 
observable by the Mariner 10 mission, in 
comparison to the dipole field. We 
submitted a paper on this study, which is 
currently in press [Stanley et al., 2005]. 
The results predict the poloidal field power 
and structure, and these are observations 
that can be tested during the 
MESSENGER mission. 
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Fig. 3. Magnetic power spectra at the surface of 
Mercury for different numerical dynamo models. 
Average power over a magnetic diffusion time vs. 
spherical harmonic degree is shown in (a) and vs. 
spherical harmonic order is shown in (b). Models 
with different inner to outer core radius ratios (r i0 ) 
and modified Rayleigh numbers (Ra m ) are shown: r i0 
=0.35, Ra m =18000 (red stars), r i0 =0.8, Ra m =25000 
(black circles), r jo =0.8, Ra m =40000 (blue squares) 
and r jo =0.9, Ra m =60000 (green diamonds). 
Differences can be seen between the models: The 
thicker, Earth-like shell thickness model (red stars) 
contains less power in degrees 3 and higher than 
thinner models, the two thin shell models with 
convection occurring both inside and outside the 
tangent cylinder (blue squares, green diamonds) 
have higher degree 3 components than degree 2 
components unlike the other models, and the thin 
shell model with convection occurring only outside 
the tangent cylinder (black circles) appears to have 
less power in axisymmetric modes (order 0) than 
non-axisymmetric modes. For more information on 
these models, see Stanley et al. [2005] (Appendix 
1 )- 
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BUDGET SUMMARY for year 3 

For period from 4/1/2004 to 03/31/2005 

• Provide a complete Budget Summary for year one and separate estimated for each subsequent year. 

• Enter the proposed estimated costs in Column A (Columns B & C for NASA use only). 

• Provide as attachments detailed computations of all estimates in each cost category with narratives as 
required to fully explain each proposed cost. See Instructions For Budget Summary on following page for 
details. 


1 . Direct Labor (salaries, wages, and 
fringe benefits) 

2. Other Direct Costs : 

a. Subcontracts 

b. Consultants 

c. Equipment 


d. Supplies 1,103 

e. Travel 2,205 

f. Other 22,628 

3. Facilities and Administrative Costs 45,783 

4. Other Applicable Costs : 

5. SUBTOTAL-Estimated Costs 135,000 


6. Less Proposed Cost Sharing (if any) 

7. Carryover Funds (if any) 

a. Anticipated amount : 

b. Amount used to reduce budget 

8. Total Estimated Costs 135,000 XXXXXXX 


| NASA USE ONLY | 
ABC 

63,281 


9. 


APPROVED BUDGET 


XXXXXX XXXXXXX 


